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Abstract 

This is the manual of the specialized event generator TopReX 3.25. The gen- 
erator provides the simulation of several important processes in pp and pp col- 
lisions, not implemented in PYTHIA (yet). Some of these processes include 
t-quarks whose spin polarizations are taken into account in the subsequent de- 
cay of the t-quarks. Several non-SM top quark decay channels are included, 
too. All calculated subprocesses can be accessed from PYTHIA as external 
processes. In addition, TopReX can be used as stand alone event generator, 
providing partonic final states before showering. In this mode the control of 
the event generation is taken by TopReX itself. A few simple examples of 
main routines, which show how to use TopReX in the different modes are 
discussed. 
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Introduction 

The event generator TopReX in its present version (3.25) provides the generation of important heavy 
particle production processes in pp and pp collisions (t-quarks, charged Higgs, etc.), which are not imple- 
mented in the high energy physics event generator PYTHIA [|l]] (yet). The subprocesses implemented 
in the present version of TopReX are listed in Table 
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Table 1: TopReX processes and corresponding process numbers. The user can choose a specific process by passing, among 
others, the wanted process number IPROC to the initialization subroutine TOPREX (FRAME, BEAM, TARGET, IPROC, 
ECM) . 

The polarization of final t-quarks is calculated and taken into account in the subsequent decay 
of the t-quarks. For processes with r-leptons in the partonic final state (originating from top decays or 
directly produced by W or bosons) the polarization of the r-leptons is calculated, too. The decay of 
the polarized r-leptons is treated by the TAUOLA package [^. Thus TopReX provides the generation 
of hard subprocesses as listed in Table |I| with subsequent decays of heavy particles, given in Table ^. 

t ^ bW+, ^ bH+, ^ qj, qg, qZ 
W,Z ^ //', f = q,l,v 
^ //', f = q,l,y 

r" Vrl'V, Vr1T-,...TAl}0LA 

Table 2: Decay channels of heavy particles included in TopReX . 

After the generation of a hard subprocess, TopReX returns for any generated event the following infor- 
mation: 

• the differential cross section value, 

• flavours and momenta of the initial and final state partons, 

• colour flow information, 

• parton shower arrangements (shower pairs/partners and scales), 

which is stored in the PYTHIA common block /PYUPPR/. Subsequently, PYTHIA can be used for the 
fragmentation of quarks and gluons into jets, followed by the hadronization and the decay of resonances. 
Initial and final state radiation as well as multiple interaction models can be used as usual. In analogy 
to PYTHIA internal processes all TopReX subprocesses can be accessed by PYTHIA through the call 
of an external process (one at a time). After the simulation of a hard subprocess in TopReX through a 
PYTHIA call of the PYUPEV subroutine, the standard PYTHIA common block /PYJETS/ is filled with 
all necessary information about the flavour and momenta of the initial and final state partons, the colour 
flow and parton shower arrangements and scales. 
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There are two alternative ways in which TopReX can be used. Fully integrated into PYTHIA 
by linkage of the TopReX package as matrix element library, respecting the programming structure of 
PYTHIA in such a way that the impact for the user has been minimized or as stand-alone program 
avoiding the overhead of PYTHIA if it is not needed. While the former is more practical for the usage 
of already existing analysis routines and interfaced detector simulations, the latter provides a complete 
analysis environment with the possibility of booking and filling histograms. 

Before the usage of TopReX will be discussed in more detail a closer look at the data structure and 
program flow gives useful insights in the functionality of the program. 

1. General structure of TopReX package 

The structure of TopReX and the corresponding program flow is shown in fig. [ij 

If needed, the user can modify PYTHIA parameters first. This can be done in the main . f routine 
or in a separate user's subroutine (e.g. subroutine TXPYIN, see appendix 2). Two examples of main . f 
routines and a corresponding subroutine to set PYTHIA parameters are presented in the appendices 1 
and 2. 

Then the subroutine TOPREX (FRAME, BEAM, TARGET, IPROC, E CM) has to be called, 
in which all variables are identical to those used in the call of the PYINIT subroutine except the integer 
variable IPROC, which specifies the number of the TopReX process. The character variable FRAME is 
foreseen to pass the information of the rest frame in which the generated event has to be evaluated. The 
colliding particles have to be specified in the character variables BEAM and TARGET and the center of 
mass energy of the collider has to be given in the variable ECM. These are all parameters, needed to run 
TopReX . Only in the case of very specific processes with Flavour Changing Neutral Current (FCNC) 
interactions the appropriate parameters should be specified in the subroutine FCNC_INI. 

After the initialization of TopReX the subroutine PYINIT has to be called for the initialization of 
PYTHIA . Subsequently events can be generated. This is typically done in an event loop as shown in the 
example program main . f. At this place the user may call its own analysis routines (see appendix 2). 

PYTHIA can access the processes of TopReX as external processes through a call to the subrou- 
tine PYUPEV. This routine, in turn, calls the appropriate TopReX routines, needed to evaluate matrix 
elements squared, Lorentz invariant phase spaces (LIPS), parton luminosities (quark and gluon distribu- 
tions of the beam particles) as well as to provide colour flow information and to set up parton shower 
arrangements and scales of the hard subprocess. Finally, all information of the generated event is stored 
in PYTHIA common block /PYUPPR/. Returned from the PYUPEV subroutine, PYTHIA apphes its 
standard procedure of jet fragmentation, hadronization and resonance decays. The simulated events are 
written to the standard common block /PYJETS/. 

As default, control printout is sent to the terminal. In addition a second possibility to store the 
printout in a separate file (toprex . info) is provided. In this case, the user should call the subroutine 
TXWRIT (-6) in the main . f routine before any initialization. The parameter value '-6' means, that 
the file toprex . info will be assigned to the Logical Unit Number (LUN) 6. 

As described above, all parameters can also be specified in two sepai^ate routines (main . f or 
TXPYIN and TopReX). In addition, the user can use a special file run . dat to read in input parameters 
(see appendix 3 for details). 

2. TopReX linkage to PYTHIA 

Along with the TopReX package comes a GNUmakefile which links the necessary libraries and the 
source code together and produces an executable. The program package was tested under Solaris, Dig- 
ital Unix, and Linux and should also run on other platforms. The libraries PYTHIA and cernlibare 
assumed to be accessible under the path /cern/pro/lib. If this is not the case the user has to change 
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the library search path accordingly. The TAUOLA library is provided as source code and will be com- 
piled as separate library. The TopReX source code, the slightly modified PYTHIA subroutines PYDECY, 
PYEVNT, PYSHOW, and PYUP IN (for PYTHIA version 6.136 and newer) are compiled together with the 
example main program listed below. To take into account the different cases of the different versions of 
PYTHIA , two targets which produce executables are provided in the GNUmakefile. For PYTHIA before 
version 6.157 the executable has to be produced by the command gmake toprex . exe while gmake 
TX . exe has to be used otherwise. If the user is already familiar with the event generator PYTHIA , the 
TopReX package may easily be integrated in existing program structures (user main program, analysis 
routines and detector simulation). 

The user has to keep in mind to compile and link all routines of the TopReX package beside his 
routines. To achieve this, either the user program file names have to be added in the provided GNUmake- 
file, e.g. they may be appended to the USEROB J variable or all TopReX related file and library names 
have to be included in the user makefile. 

3. Interface and program flow 

Global TopReX parameters can be accessed and changed via the common block TXPAR: 



COMMON/TXPAR/ Ipar(200), Rpar(200) 



Integer values like the TopReX process number to be chosen by the user are stored in the array I par 
as indicated in Table |]. Floating point values like the hadronic centre-of-mass energy available for 
collisions are stored in the array Rpar. A complete list of the parameters and their default values is 
given in appendix 4. As a minimal requirement the number of a subprocess has to be specified. Further 
parameters like the hadronic centre-of-mass energy, particle masses etc. are taken from PYTHIA . Two 
examples how to use TopReX in conjunction with PYTHIA are given in appendices 1 and 2. 

In analogy to the PYTHIA common block /PYDAT3/, which offers the possibility to switch 
on/off specific decay channels of particles after the generation of the hard subprocess, the TopReX 
common block /TXRDEC/ provides this possibility for W^, H^, and i-quark decays which are part of 
the hard subprocess. For convenience, the numbering scheme is adopted to that of PYTHIA. Numbers 
of the Individual Decay Channels (IDC's) for W and Z-bosons are identical to those of PYTHIA (see 
PYTHIA Manual IDC numbers for are identical to those of the boson with one additional 
channel (IDC=21): bt* bWb. For top quark decays the IDC's are defined as follows: 

1. (IDC = 1) SM decay, t bW 

2. (IDC = 2)t ^ bW decay via SM and beyond SM interactions 

3. (IDC = 3)t ^ bH^ 

4. (IDC = 4,. .,11) correspond to top quark decay via Flavour Changing Neutral Current interactions 
The common block /TXRDEC/ is defined as follows: 



COMMON/TXRDEC/ MID (4, : 30) , BRF(4,30), FID(4,30,5), BRS(4,5) 

with the integer parameter array MID and the double precision parameter arrays BRF, FID and BRS, 
explained below. 

Purpose: to access particle decay data and parameters. The first index I (1=1,2,3,4) of the arrays 

corresponds to t-quark, W, H, and Z-boson decays, respectively. 
MID (I , IDC ) : on/off switch for IDC (meaning identical to PYTHIA parameter MDME ) 
MID (I, 0) : total number of allowed IDC 
= - 1 coupling switched off 
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= channel is switched off for decay, but it contributes to the total width 

= 1 channel is switched on 

= 2 channel is switched on for particle but off for anti-particle 

= 3 channel is switched on for anti-particle but off for particle 
= 4,5 are not used in current version 3.25 of TopReX . 
BRF { I , IDC ) : branching fraction for given IDC 

FID{I,IDC,j) : contains the KF code for decay products of a given IDC, 

meaning identical to PYTHIA parameter KFDP 
BRS { I , IDC ) : sum of branching fractions for different decay groups 

4. Description of hard processes 

Some comments on the hard processes of TopReX , listed in table |l|, are made in the following. The kine- 
matics of a hard process can be specified by setting the values of the PYTHIA parameters CKIN { 1 : 4 ) 
in the main . f or TXPYIN routine. The values of the CKIN parameters are copied to the TopReX pa- 
rameters Rpar (101 : 104) in the subroutine TOPREX. Three processes (Ipar { 1 ) =1, 7, 10) 
are singular at p_\_ — > 0. To avoid the singularities a conventional kinematical cut on (px) has to be used 
(a minimal value of 20 GeV is recommended): 

CKIN { 3 ) = 2 . dO (copied during initialization to the TopReX parameter Rpar (10 3)) 

The 1^-gluonfusion processes (Ipar (1 ) = 1 , 2), W^* ^ rz^ (Ipar ( 1 ) = 6,7)andff^* — 
TV (Ipar ( 1 ) = 9,10) need some comments concerning the problem of double counting, given be- 
low. Following the work of [^] a simple method for the generation of such reactions is applied (see also 

o Single top Production: 14^-gluon fusion, Ipar ( 1 ) =1, 2 

Two different subprocesses can be used for the generation of single top production in 1^-gluon fu- 
sion processes: 2 ^ 2 (Ipar (1 ) = 1) with CKIN (3) = . and 2 ^ 3 (Ipar (1) = 2) with 
CKIN ( 3 ) = 2 0.0. The additional &-quark in the first kind of events appears in the initial state and can 
be found in the PY JETS list. Such an event will be accepted for further analysis, if the transverse momen- 
tum of this additional 6-quark does not exceed some threshold p_lo(^) (typically of the order ~ 10 GeV). 
An event of the second kind {qg — > q'tb process) will be accepted if the transverse momentum of the 
final 6-quark from the hard process will be greater than the p_lo(&) value. 

oW* ^ TV {IpsLT {!) = 6 ,1) and H^* ^ TV {Ipar {!) = 9 , 1 0) processes: 
These processes are generated in the same way as for the single top production process through the 
1/l^-gluon fusion mechanism. However, the user should examine the transverse momentum of the inter- 
mediate W* or i7* -boson, here. For W* production processes it is recommended to choose threshold 
values above 

P±o{W) « 40 GeV 

For the charged Higgs production the p±o{H) value depends on the mass of the charged Higgs and 
should be figured out by the user. 

For processes with a charged Higgs (H^) the default values for M{H^) and tan/3 are given 

below: 

M(i7=t) = 300 GeV Rpar (51) 
tan/3 = 50 Rpar (55) 

To change these values the corresponding parameters Rpar ( 51 ) , Rpar ( 55 ) have to be modified in 
the subroutine TOPREX. 

o Top anti-top quark pair production: tt —>■ hW^hW~ —>■ 6 fermions 

The processes Ipar(l)=2 0,21,22 provide the production in the Breit-Wigner approach with off- 
shell W bosons [Q] (see for details about the top quark production). The tt processes (Ipar(l)=50, 



4 



51, 52) in the on-shell approximation [§] provide the possibihty to switch on an intermediate neutral 
Higgs boson which becomes resonant in the case of a Higgs mass above the tt production threshold. 
Below, the ti production via a Higgs resonance is switched off. The matrix elements calculate the co- 
herent contribution of the resonant Higgs boson and the non resonant background. They are also valid 
in supersymmetric extensions of the standard model in case of squark masses above 400 GeV, since the 
contributions of such heavy squarks are negligible [Q]. The scalar and pseudoscalar Yukawa couplings 
to top quarks are described by the parameters Rpar (81) = a and Rpar (82) = a respectively. The 
standard model couplings of the Higgs boson to the W and Z vector bosons are realized with a multi- 
plicative factor, the relative coupling strength Rpar (80) = gvv^ which depends as the reduced Yukawa 
couplings in a two Higgs doublet model on the ratio of the vacuum expectation values of the Higgs dou- 
blets. In the ti matrix elements used here this factor is important, to take into account the total decay 
width of the Higgs boson, which is also printed as control output after initialization of the processes. The 
user has to provide the values for the Higgs boson couplings. Several examples are listed in Table |3[ 

Rpar (80) Rpar (81) Rpar (82) 
Standard model Higgs boson: gw = a = 1, a = 

Pseudoscalar Higgs boson: gw = 0, a = 0, a = 1 

CP violating Higgs boson: a / 0, a / 

Table 3: Examples of the Higgs boson couplings Rpar ( 8 ) = Qvv, Rpar ( 81 ) = a and Rpar ( 82 ) = a of the tt production 
processes Ipar(l)=50,51,52. 



5. Comparison of parton distributions 

To demonstrate the reliability of the processes implemented in TopReX , they are compared to the default 
processes of PYTHIA 6. 1. for pp collisions at ^/s = 14 TeV (LHC). 

In the upper plot of figure |2| the distribution of the top quark mass is shown for the two event 
generators. The 2 — > 6 matrix elements in the Breit-Wigner approach are in good agreement with the 
2 — > 2 tt matrix elements of PYTHIA 6. 1. The lower plot shows that the p± spectrum of the top quarks 
are in agreement for the different event matrix elements of the compared generators. 

In figure ^ the spin correlation between the top quark decay products is shown for the tt produc- 
tion in the di-leptonic decay channel. The angles of the two leptons are evaluated in the helicity basis 
convenient at the LHC [jl4|], [p^]. While the 2^2 matrix elements of PYTHIA cannot include the spin 
correlation, the 2—^6 matrix elements of TopReX do. Conclusively, the asymmetry coefficient of the 
spin correlation obtained using default PYTHIA gives vanishing values. In contrast the standard model 
prediction, reproduced by TopReX to leading order, yields the asymmetry coefficient Asm = 0.32 using 
the parton densities of CTEQ4L. 

Figure ^ shows a comparison of the tt spin correlation with help of the helicity basis for the matrix 
elements in the Breit-Wigner (left) and on-shell (right) approach, implemented in TopReX . The contri- 
bution to this correlation of the gluon gluon fusion (top) and the quark anti-quark annihilation processes 
(center) is given separately. The predictions of the matrix elements in the Breit-Wigner approach and in 
the on-shell approach agree quite well. 

6. Conclusions 

The specialised event generator TopReX , which provides the simulation of several heavy quark produc- 
tion processes in hadronic collisions has been described. As cross check many comparisons of parton 
distributions between the matrix elements implemented in PYTHIA and TopReX have been applied. The 
agreement is very good except in cases Uke the spin correlation where default PYTHIA is not able to 
reproduce the standard model predictions. 



5 



Acknowledgements 

Great thanks to S. AbduUin, W. Bernreuther, E. Boos, K. Datsko, D. Denegri, V. DroUinger, G. Fliigge, 
U. Heintz, V. Ilyin, R. Mehdiyev, A. Miagkov, D. Murashev, A. Nikitenko, V. Obraztsov, R. Schulte, 
S. Stepanov, and O. Yushchenko for fruitful discussions and support during the development of TopReX 
and this manual. 



6 



References 

[1] T. Sjostrand and M. Bengtsson, Comput. Phys. Comm. 43, 367 (1987); T. Sjostrand, Comput. Phys. 
Comm. 82, 74(1994) 

[2] S. Jadach, Z. Was, R. Decker and J. H. Kuhn, Tauola 2.4 Comput. Phys. Commun. 76 (1993) 361. 

[3] S. Abdullin, A. Khanov, N. Stepanov, CMS TN/94-180 (1994). 

[4] R. Kleiss, W. J. Stirling, Z. Phys. C 40, (1988) 419. 

[5] M. Beneke et al., in [0] 419-529 [ |hep-ph/0003033| ] . 

[6] W. Bemreuther, M. Flesch, P Haberl, Phys. Rev. D58, (1998) 114031. 

[7] S. Dawson, A. Djouadi, M. Spira, Phys. Rev. Lett. 11, 16 (1996). 

[8] E.E. Boos, L.V. Dudko, V.I. Savrin, CMS NOTE-2000/065 (2001). 

[9] S. Slabospitsky, CMS-Note (2002) in preparation. 
[10] Y. P Gouz and S. R. Slabospitsky, Phys. Lett. B457 (1999) 177 [ ]hep-ph/981 1330| l. 
[11] S. R. Slabospitsky and G. Azuelos, Phys. Lett. B480 (2000) 105 [ |hep-ph/0003024| |. 
[12] L. Sonnenschein, CMS-Note 2001/001 (2001). 
[13] R. Brun and D. Lienart, Version 4," CERN-Y250. 
[14] G. AltarelU, M. L. Mangano et al, CERN-2000-004. 

[15] W. Bemreuther, A. Brandenburg, Z. G. Si, P Uwer, DESY-01-085, PITHA 01/06, TTPOl-16. 



7 



Appendix 1. 

The example main program shows the usage of a particular (PYTHIA external) TopReX subprocess. 
PYTHIA parameters like particle masses, etc. and the TopReX process number (IPROC parameter) are 
specified in the MAIN_1 routine shown below. To access any TopReX subprocess the parameter IPROC 
has to be set to the desired subprocess (here IPROC = 50 specifies the tt production in the on-shell 
approach). The choice has to be made before the call of the TOPREX subroutine. The further initialisation 
of PYTHIA is done by a call to the subroutine PYINIT. 

PROGRAM MAIN_1 
IMPLICIT NONE 

CHARACTER* 8 FRAME, BEAM, TARGET 

DOUBLE PRECISION ECM ! CMS energy 

INTEGER IPROC ! number of TopReX process 

INTEGER I 

EXTERNAL PYDATA 

INTEGER MSTU, MSTJ,KCHG, MDCY,MDME, KFDP, MSEL, MSELPD, MSUB 
INTEGER KFIN, MSTP, MSTI 

DOUBLE PRECISION PARU, PARJ, PMAS, PARE, VCKM, BRAT, CKIN 

DOUBLE PRECISION PARP, PARI 

COMMON /PYDATl/ MSTU ( 2 ) , PARU ( 2 ) , MSTU ( 2 ) , PARU ( 2 ) 
COMMON /PYDAT2/ KCHG ( 5 , 4 ) , PMAS ( 5 , 4 ) , PARE ( 2 ) , VCKM ( 4 , 4 ) 
COMMON /PYDAT3/ MDCY ( 50 , 3 ) , MDME ( 4 , 2 ) , BRAT ( 4 ) , KFDP ( 4 , 5 ) 
COMMON /PYSUBS/ MSEL, MSELPD, MSUB ( 500 ) , KFIN ( 2 , -4 : 4 ) , CKIN ( 2 00 ) 
COMMON /PYPARS/ MSTP ( 2 00 ) , PARP ( 2 00 ) , MSTI ( 2 00 ) , PARI ( 2 ) 



FRAME 
BEAM 
TARGET 
ECM 

IPROC = 



' CMS' 
' P' 
'P' 

14000 .dO 



50 



! LHC cms energy (in GeV) 
! tt" production 



PMAS (6,1) = 
PMAS (5,1) 

PMAS (23,1) = 

PMAS (24,1) = 

PMAS (25,1) = 

MSTP (51) = 7 



175. dO 

4 . 8d0 
91 . 187d0 
80 . 41dO 
400 .dO 



top mass 
bottom mass 
ZO mass 
W mass 
Higgs mass 



! PDF: CTEQ5L 



CALL TOPREX (FRAME, BEAM, TARGET, IPROC, ECM) 
CALL PYINIT (FRAME, BEAM, TARGET, ECM) 



DO I = 1,10 

PRINT*, ' Event number 1=' , I 

CALL PYEVNT ( ) 

CALL PYHEPC(l) 

IF (I.LE.2) CALL PYLIST(2) 

at this place analysis routines and/or 
detector simulations may be called 



ENDDO 



END 
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Appendix 2. 

Contrary to the previously described MAIN_1 routine, this example considers the possibility to set all 
PYTHIA parameters in a special TXPYIN subroutine. In addition, a few parameters are read in by the 
TXRINT subroutine from an external input data file, named run . dat (see appendix 3 for details). An 
example of a special routine for user's analysis (USRPRO) is described, too. 

PROGRAM MAIN_2 

IMPLICIT NONE 
DOUBLE PRECISION Ecm 
INTEGER IPROC 

INTEGER Ntot, N, MODE, lER 

Ecm = 14000. dO ! CMS energy (in GeV) for LHC option 

call TXWRIT(-6) ! LUN=6, open file with output information 

call TXRINT (' iproc=' , iproc, ier) ! no. of TopRex process 
call TXRINT (' Ntot=' , Ntot, ier) ! number of events to be generated 

call TXPYIN (IPROC, Ecm) ! all PYTHIA will be specified there 

* for compilation user' s program WITHOUT TopRex package comment next line 
if (Iproc. ge. 1) call TOPREX (' CMS ',' P ',' P ' , IPROC, Ecm) ! TopRex init . 

call PYINIT ('CMS' , 'P' , 'P' , Ecm) 

call PYSTAT(4) 

c... call analysis/detector simulation routine for initialisation 

mode = -1 

call USRPRO (mode, Ntot) 
end of initialisation 
DO N = 1, Ntot 

if (mod (N, 1000) . eq. 1) write(*,*)' event no=' , N 

call PYEVNTO 

call PYHEPC (1) 

c... call analysis/detector simulation routine for running 
mode = 

call USRPRO (mode, N) 
ENDDO 

call PYSTAT(l) ! Brief statistics output from PYTHIA 

*... call analysis/detector simulation routine for closing 
mode = 1 

call USRPRO (mode, N) 

STOP 
END 

The following example routine TXPYIN is used to set PYTHIA parameters. 

SUBROUTINE TXPYIN ( IPROC, Ecm) 
IMPLICIT NONE 

*... Standard PYTHIA ( v. >= 6.1) commons for initialization. 
EXTERNAL PYDATA 

INTEGER MSTU, MSTJ, KCHG, MDCY, MDME, KFDP, MSEL, MSUB, KFIN, 
& MSTP, MSTI, MSELPD 

DOUBLE PRECISION 
& PARU, PARJ, PMAS, PARE, VCKM, BRAT, CKIN, PARP, PARI 
COMMON /PYDATl/ MSTU ( 2 ) , PARU ( 2 ) , MST J ( 2 ) , PARJ (2 ) 



CMS energy 

no. of Toprex process 

internal variables 
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COMMON /PYDAT2/ KCHG ( 5 , 4 ) , PMAS ( 5 , 4 ) , PARF ( 2 ) , VCKM ( 4 , 4 ) 

COMMON /PYDAT3/ MDCY (500, 3) , MDME (4000,2), BRAT (4000) , KFDP (4000,5) 

COMMON /PYSUBS/ MSEL, MSELPD, MSUB ( 500 ) , KFIN ( 2 , -4 : 4 ) , CKIN ( 2 GO ) 

COMMON /PYPARS/ MSTP ( 2 00 ) , PARP ( 2 00 ) , MSTI ( 2 00 ) , PARI ( 2 ) 



INTEGER IPROC, I 



MSTP (51) = 8 

masses . . . , 

PMAS(6,I) = 175. dO 

PMAS (5,1) = 4.8d0 

PMAS(24,1) = 80.4IdO 

PMAS(25,1) = 400. dO 

.KINEMATICS 



PDF CTEQ5M1 (PYTHIA6 . 1 . . ) 

top mass 
bottom mass 
W mass 
Higgs mass 



CKIN (3) = 2 0.d0 ! PT HAT LOW CUT 

CKIN (4) = 200.d0 ! PT HAT UP CUT 



if (IPROC .NE . 0) return ! IPROC > 0, TopRex process 
standard PYTHIA process is chosen, no TopRex call ! 



MSEL = 
MSUB (16) =1 



! ffbar — > gW+ or gW- 



RETURN 
END 



This example routine is used for further fast simulation of CMS detector response (by call of CMS JET routine) and 
for writing of NTUPLE (by call of TX_NT routine). 

SUBROUTINE USRPRO(MODE, NV) 

IMPLICIT NONE 

INTEGER MODE, NV ! internal variables 



if (MODE.EQ.-l) THEN 
initialisation 

call CMSJET (MODE) 
elseif (MODE .EQ . 0) THEN 
call CMSJET (MODE) 

elseif (MODE .EQ . 1) then 

call CMSJET (MODE) 
ENDIF 
END 



! reading and initialisation of CMSJET 

! running 

! running of CMSJET 

! closing files, etc 
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Appendix 3. 



Tools to read in parameters 

The input data file run . dat has a simple format. The data card will be ignored, if it is: 

• a totally blank card 

• a card with the first non-blank character (a comment card) 

• a card without '=' character 

The first item in the card with data is a key of a variable. The '=' character should be expUcitly 

given behind this key. Then one or several numerical values have to be given. The rest of the card is 
considered as comment. Two examples of such cards are given below. 

• comment 

iproc =21 ! no. of TopRex process (21 = gg -> tt" production) 

filout = myout.dat ! name of the output data file 

• New Physics scale (1 TeV) 
Rpar(60) = 1000.0 ! in GeV 

In this example the items 'iproc =', 'filout =' and 'Rpar (60) =' are the keys, while 
'21', 'myout.dat' and '1000.0' are the corresponding values to be read. Four different subroutines are 
provided to read these data: 



call TXRCHR(fkey, 'filout', fname, lER) 

call TXRINT(fkey, 'Iproc =' , iproc, lER) 

call TXRFL8 (fkey, 'rpar (60) =' , Rpar (60), lER) 

call TXRFL4 (fkey, 'name =' , realvalue, lER) 



character value 
integer value 

double precision value 
single precision value 



where the second parameter is a key-string, the third parameter is the value of the key variable, and 
lER is an error flag (lER = 1 means that the subroutine is not able to find appropriate data in the file 
' run . dat ' ). 

The first parameter, fkey, is the key for the file, from which the parameters are to be read (the 
' run . dat ' file in this example): 

character *6 fkey ! key for file with input data 
data f key/ ' rundat ' / 

The user can also read several numerical values from one card. As an example, we present the 
card with the key of 'H+- boson', with two numerical values, '220.' and '25.' being the mass of the 
i?^ -boson and the tan P parameter. 

* mass and tan (beta) parameters for H+- boson process 
H+- boson = 220 . 25 . 

To read this card a character string (here STROUT) should be described and the TXRSTR subrou- 
tine is used to read these two values: 

CHARACTER *80 STROUT ! string to be used for data input 



call TXRSTR('h+- boson =' , STROUT, lER) 

read ( STROUT, *) Rpar (51), Rpar (55) ! mass and tan (beta) 



For string-key values in the run . dat file and in the parameters of the subroutines TXR . . . 
lowercase and uppercase letters are allowed (the read in routines are not case sensitive). Blank characters 
are allowed, too. 
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Appendix 4. 

The parameters of the /TXPAR/ common block are explained in the following. The default values of 
the parameters are given in the parentheses (D = ...) 



COMMON/TXPAR/ Ipar(200), Rpar(200) 



Purpose: 1 


Ipar 


(1) : 


Ipar 


(2) : 


Ipar 


(4) : 


Ipar 


(5) : 


Ipar 


(6) : 


Ipar 


(7) : 


Ipar 


(8) : 


Rpar 


(1) : 


Rpar 


(2) : 


Rpar 


(3) : 


Rpar 


(4) : 


Rpar 


(10) 


Rpar 


(11) 


Rpar 


(12) 


Rpar 


(13:1 


Rpar 


(18:1 


Rpar 


(20 ) 


Rpar 


(21:2 


Rpar 


(30) 


Rpar 


(31) 


Rpar 


(32) 


Rpar 


(41:4 


Rpar 


(45 : 4 


Rpar 


(51 :5 


Rpar 


(60) 


Rpar 


(61 : 7 


Rpar 


(81) 


Rpar 


(82) 


Rpar 


(83) 


Rpar 


(101) 


Rpar 


(102) 


Rpar 


(103) 


Rpar 


(104) 



.am in GeV. 
^2^inGeV2. 



TopReX 's process number 
number of events to be generated 
continuation flag 

number of entries used for estimation of SigMax (D = 200000) 
number of entries used for estimation of |Mp (D = 100000) 
KF code of Q-quark for WQQ process (D = 5, (5=6-quark) 
quark masses, (D=l : RPP values), 2 : PYTHIA values 

^/S = Er 

s 

evolution scale for PDF evaluation 

evolution scale for evaluation of ag 

oiQED = (1/128) electromagnetic coupUng 

e = ^/^koqed electric charge 

Fermi constant G p 
' ) : sin cos i^w^ sin^ i^w, cos '&w^ sin 2'&w 
)) : g = e/ sin-dw, = e/ sm2'dw 

0.38939 transformation l/GeV^ to mb 
i) : TT, 2^, (27r)3, (27r)4, (27r)6 

SM top quark decay width (t bW) 

total top quark decay width (including non-SM interactions) 

partial top quark decay width for t — > bH^ 
1 ) : W-boson parameters, Mw, M^, r{W), Mw ■ ^{W) 
] ) : Z-boson parameters, Mz, Ml, T{Z), Mz ■ T{Z) 
]) : boson parameters, Mh, Mjj, V{H), MT, tan /3, cot /3 

(D=1000.dO) New Physics scale, A = 1 TeV 
) ) : anomalous FCNC top quark couplings (FCNC_INI routine) 

gvv^ relative coupling of Higgs boson to W and Z vector bosons 

(processes Ipar(l) = 50, 51, 52 only) 

a, scalar Yukawa coupling of Higgs boson to top quarks 

(processes Ipar(l) = 50, 51, 52 only) 

a, pseudoscalar Yukawa coupling of Higgs boson to top quarks 
(processes Ipar(l) = 50, 51, 52 only) 
\/smin for hard process (D = 20 GeV), identical to CKIN ( 1 ) 
y/smax for hard process (D = ^/s), identical to CKIN ( 2 ) 
minimal value for p±(D = 5 GeV), identical to CKIN ( 3 ) 
minimal value for p± (D = \fs/2), identical to CKIN ( 4 ) 



12 



Program flow chart 



Program 



MAIN 



Initialization 



I PYTHIA 
I parameters 



Event loop 



analysis 
routines 



End 



PYTHIA 



PYINIT 



PYEVNT 



PYRAND 



common /PYJETS/ 



TopReX 



TOPREX 



TX_TIT 



TXWEAK 



PYUPEV 



\M\'^0LIPS 



\ 

common /PYUPPR/ 



Fig. 1 : Program flow chart. The main program, stored in the physical file main.f, initializes TopReX and PYTHIA . The initial- 
ization of TopReX consists predominantly of the passage of electroweak parameters from PYTHIA to TopReX followed by the 
estimation of the maximal value of the differential cross section for the chosen hard scattering process. After the initialization 
the generation of scattering events takes place in the event loop of the main program through a call to the PYTHIA subroutine 
PYEVNT. The subroutine PYUPEV (whose PYTHIA dummy version is overwritten by TopReX ) is called subsequently. 
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Fig. 2: The top quark mass (upper plot). The distribution of the 2 6 ti matrix elements of TopReX in the Breit- 
Wigner approach shows agreement with the prediction of the 2-^2 matrix elements of PYTHIA 6.1. The natural 
width of the top quark mass is fitted with a Breit-Wigner function. The lower plot shows the pi_ spectrum of the top 
quarks. The prediction of PYTHIA 6. 1 coincides with the on-shell approach of TopReX . The Breit-Wigner approach, 
used here with a completely different scale, gives quite similar results. 
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PYTHIA 6.1 
partonic final state 
87%gg+13%qq 



LHC 14TeV 
CTEQ4L 



•t—i 
c 

> 

0) 




.-1 



-1 



PYTHIA 6.1 + M.E. Slabospitsky 
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Fig. 3: tt spin correlation in the helicity basis. The prediction of the 2^2 matrix elements of PYTHIA 6.1 show 
no correlation between the two leptons (upper plot). In contrast, the standard model predicts the correlation shown in 
the lower plot. This result is obtained by the 2 — > 6 matrix elements implemented in TopReX . In the Breit-Wigner 
approach (shown here) the asymmetry coefficient reads Asf^ = 0.32. In the on-shell approach the coefficient amounts 
to = 0.33. 
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Fig. 4: ti spin correlation of TopReX matrix elements in the helicity basis. To the left the prediction of the 2^6 
tt matrix elements in the Breit-Wigner approach is shown. To the right, the same is shown for the 2 ^ 6 tt matrix 
elements in the on-shell approach. At the top only the contribution of the gluon gluon fusion is plotted. In the center 
only the contribution of the quark anti-quark annihilation is given and at the bottom both production processes are 
taken into account. The corresponding asymmetry coefficients are indicated in the plots. The Breit-Wigner and on- 
shell approaches are in very good agreement. 
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